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ABSTRACT 

The subject of this paper is the mechanical behavior of rectangular plates subjected to 
a combination of axial compression and lateral pressure. Displacements and stresses are 
obtained from a Fortran code based on the von Karman plate equations. The effects of 
various boundary conditions, nonlinearities, and imperfections are included. 
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INTRODUCTION 



Fatigue life predictions require knowledge of the stresses in a ship under its operating 
conditions [Sikora, Dinsenbacher, and Beach (1983)]. Stiffened plates are a basic structural 
component of ships and submarines. The mathematical equations governing the deformation 
of thin elastic plates and methods for solution of these equations are well-known [Timoshenko 
and Woinowsky-Krieger (1959), Szilard (1974), and Hughes (1983)]. The objective of our 
work is to use these known mathematical methods to predict the stresses in typical ship 
plating. 

Plates are subjected to axial tension and compression from longitudinal bending of a ship 
due to wave loads. In addition, plates on the bottom are subjected to lateral water pressure/ 
Deck plating boundary conditions may be taken to be simply-supported , whereas bottom 
plating boundary conditions are more closely approximated as clamped. Deck plating may 
be considered to have an initial geometric imperfection, whereas bottom plating is bowed in 
by the water pressure. 

We consider a rectangular plate with length a, width b, and thickness t. The structure 
is subjected to axial force F, and possibly uniform lateral pressure p (force per unit lateral 
area of the plate). 
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PLATE EQUATIONS 



The basic differential equations of nonlinear shallow shell theory are 

D V ' w = P + ( w 0 + w), xx -\-^, xx (w 0 + w), yy — 2 4* ixy ( w o + u: ) -xy ( 1 ) 



— V 4 O = [(w 0 + w),l y ( Wq + w), xx (w 0 + w),yy } ~ Oo,ry ~W 0 <xx W 0 , yy ) (2) 

Here (j\y,z) are Cartesian coordinates measured from the center of the plate. 

D = Et 3 / 12(1 — i/ 2 ) is the bending stiffness, E is Young’s modulus, and v is Poisson’s ratio. 

Wq(x, y) is the initial geometric imperfection of the plate midsurface in the z-direction, while 

w(x,y) is the normal displacement. 4>(x,y) is the Airy stress function. Commas denote 

partial differentiation with respect to x or y; e.g, 1 

A d 4 w d 2 w d 4 w 

V W = W ixxxx +-U>,xxyy +W ,yyyy = + “^ 2^2 + 

We let u(x,y) and v(x,y) denote the displacements of the plate midsurface in the x and 

y directions, respectively. These in-plane displacements can be related to w and <t> by using 

the strain-displacement relations 

e x = u, x +w 0 , x w, x +^w,l 



C y = V,y+W 0 , y W, y +\w,l 

txy — +U, X +u;0,x -\-Wo,y W, x -\~W, X W,y 

and the constitutive relations 





Et . 


N x = 


1 _ u 2 ^ x + Ut v) 




Et 


N y = 


1 _ u 2^ e y + U€x ) 




Et 


II 


2(1 + u) €xy 



The membrane stress resultants (forces per unit length of side) are related to the stress 
function by 



:V — d) N — <I> N = —<$> 

x — ’ 1 V — ^ ixx ? 1 y xy ^ ixy 
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The bending stress resultants (moments per unit length of side) are related to the normal 
displacement by 



M x = —D(w. xx +i/w, yy ), M y = -D(w, yy +i/w, xx ), M xy = D{\ - u)u\ 



xy 



(3) 



The bending stresses (forces per unit area of cross section) acting on the top of sections 
parallel to the yz or xz planes, respectively are 



6A/ r 



6 M v 



(J X 'i ' v — 

f 2 y t 2 



We consider three different possible sets of boundary conditions: 
Loosely-clamped 



On 


* = ±f : 


u> = w, x = X xy = 0, u — constant, f\ Z x (± 1 2 ,y)dy = F 

2 


On 


y = ±f : 


a 

w — iL\ y — N xy = 0, v = constant, fj* N y (. r, +|, )dx = 0 


Rigidly-clamped 




On 


* = ±t : 


b 

w = w, x = X xy - 0, u - constant, f\ N x (±^y)dy = F 

2 


On 


y = ±2 '■ 


W = W,y = N X y = V = 0 


Simply-supported 




On x 


= ±! : w 


b 

= w,xx=N xy = 0, u - constant, f\ N x ( + \,y)dy = F 

2 


On y 


= ±| : w 


a 

— te, yy = iV ry = 0, v = constant, /_ 2 a Af y (x, +|)dx = 0 



(4) 



(5) 



( 6 ) 



The classical buckling load of a simply-supported plate with an integral aspect ratio 



- = 1,2,3,... is given bv — F cr , where 
6 



F ct = 



4t t 2 D 
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SIMPLIFIED EQUATIONS 



For normal displacements w which are smaller than the plate thickness t, the nonlinear 
terms in (l)-(2) are negligible and we can replace (l)-(2) with the linear approximation 



F 

N x = — , N y = constant C, N xy = 0 



( 7 ) 



D \/ A w — p + N x (wq + w), xx +N y (wo + u>), 



yy 



( 8 ) 



For the boundary conditions ( 4 ) and (6) where the lateral edges y = ±- are free to expand 
or contract N y = 0, whereas for the boundary conditions (5) where the lateral edges y — ±- 



are restrained N v = 
form 



v 



F 



y — . If p 7^ 0, equation (8) mav be written in the nondimensionalized* 

b 



4^2 p ] S { 

\y 4 W= 12(1 -V 2 ) + — {wq + w),xx H ^-(w 0 + w),yy 

t cr U 



w FtF x 

where w = — — , x = — , etc. The nondimensionalized bending stresses are then 
pb 4 b 

given by 



cr x t 2 6 M x 



1 



"(^511 )} etC. 



pb 2 pb 2 2(1 — u 2 ) 

Note that in the case wq = 0 the normal displacement w increases linearly with the pressure p, 
whereas in the case p = 0 the normal displacement w increases linearly with the imperfection 
amplitude. 

We can investigate the importance of nonlinearities by substituting the linear solution 
w to (8) into the right side of (2). The additional membrane stresses S x ,S y ,S xy are then 
determined by 



P p 

Nx = -T b tS X = — + 0,^ , Ny = tSy + C = 4> + C, N X y = t S xy — ~ ^ , X y 

b b 
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V 4 ^ = [(u-’O + 11'),% -(«’ 0 + «’).« (ti’O + U’),yy ] - (ttfo.xy -ti ’ 0 ,xx U’Oiyj, ) 



(9) 



£f 



If /> 7 ^ 0, Equation (9) may be written in the nondimensionalized form 



V 4 <^) = [(W’O + _ ( U ’0 + «’)« ( t«0 + »’)-yy] ~ (^’0.|y “ "’O-xx «’0 -y y ) 

_ 4>£i 5 

where T = — — — The nondimensionalized additional membrane stresses are then given bv 



S x t 2 pb 4 — 

— — — d) — etc 

pb 2 Et*' yy 



( 10 ) 



F 



Note that the total longitudinal stress at the top of a plate cross section is — + < 7 X + S T , etc.. 

tb '• 



SOLUTION METHOD 

These equations may be solved numerically by superimposing a Fourier series particular 
solution plus Levy-type homogeneous solutions with coefficients determined by the prescribed 
boundary conditions. I his method of solution has been recently used by Kwok, Kang, and 
Steele (1991), Bird and Steele (1991, 1992), and Kang, Wu, and Steele (1993). A thoroughly 
documented Fortran code was written to implement this solution. 

Normal Displacement and Bending Stresses 

The governing equations for an isotropic plate under axial compression and lateral pressure 
load with no imperfection are from (7) and (8) 

Z7 1 zd 

N x = J, Ny = 0 or N xy = 0 

D V 4 W = P + iV x W, XX +jVy W,yy 

with the clamped boundary conditions from (4) or (5) 

Onar = ±^: w = w, x = 0 
On y = ±| : w — te, v = 0 

The solution to the problem above is obtained in several stages: 
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1 ) 1- list a solution w j for constant pressure load on a clamped plate strip with A r x = S y = 0 
is obtained. 1 his particular solution that satisfies 

v 4 w = jj ( ll ) 

is „f t he fullovving form 

“' = ^('-(t)T (i -> 

The solution w\ has zero rotation and displacement along y = ±-. From (3) the bending 
moment M y is given by 



\ly DW,yy 



IS AD 
6 2 



(l — 3 j an( l M X = l 'My 



p 

(’oeffi< lent A in (12) is determined by substituting (12) in (11): A = — -^j-. 



he ri 



airier cosine series for (12) is given by 



w x = T, COS 



mry 



where 



4 /f uny pb 4 ( - 1 12 

I t n = J, J 0 w i cos “ — “1/ = T74 I 7 — — + 



717T 



3 D y(n7r) 3 (7Z7r) 5 y 



sin 



The particular solution W\ has zero rotation but nonzero displacement along x — + — 



2) In order to make the normal displacement zero along x = + — , we need to compute a 
correction term that will achieve this goal. This new term 1 x 2 which is the solution of a 
rectangular plate with inplane resultants N x and N y is computed in the following way: A 



solution of 



f the form w 2 = ^ ( x ) cos — - is assumed, where X n is a function of x only 



and is d< 



termined from the condition that w 2 satisfies the plate equation 



D \/ TV N X TV ) XX T Xy TV ^yy 



(13) 
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wit h zero rotation and displacement equal to —w\ along x = ± — . From the condition that 
fr> satisfies (Id), we obtain the solution 

-V u (jt) = A n cosh V\x + B n sinh V\x + C n cosh u^x + D n sinh 

where U\ and //•> are solutions of 




The constants B n and D n can be taken to be zero due to the symmetry about the //-axis, and 
the remaining constants A n and B n can be determined by imposing the boundary conditions 
along x — ±— . Then uq + u>2 is a solution for a plate under lateral pressure p and in-plane 

resultants \ r and A r y with zero displacement along all the edges and zero rotation along 

a . b 

x — ± — . but nonzero rotation along u = ±-. 

6 

3) In order to make the rotation zero along y — ±- we first compute the Fourier coefficients 
of rotation on these edges and then obtain a third solution w 3 = w 3 \ + w 3 2 with moments 
distributed along the edges corresponding to the desired rotations. A solution W31 sinusoidal 
in the r-direction is: 



mirx [ cosh fi\ y cosh ^2 y 



lu 3 1 — ^2 ^ lm (OS „ l 1 b 1 b 

modd a V cosh /i 1 ^ cosh fl 2 5, 



where the parameters /i r and \ii are the roots of the polynomial 



/m7r \ 2 



° e - 1 . 



r . m 7 r\ 2 ~ 

+ A x ( — — J — N y i-i~ — 0 



This solution is obtained in a similar manner as for solution W2, except that in this rase 

mux 



IV. 



31 is assumed to be of the form te 3I = ^ Y m (y) cos and Y m {y) is determined by 



imposing the zero displacement boundary condition along y = ±-. As defined above, w 3 \ 



has zero displacement on all the edges and zero moment M y along x = ±— . The rotation 
along the edge y = - is given by 



roty = ~{w 3 i, y )| y= ! = ^ C Xm cos 

TTl odd 

^ ^ lm COS 



mnx 

a 

mnx 



f . b , />' 

— /ii taidi - + /i 2 tann /r 2 - 






(15) 



lm 



The rotation along the edge x = — is 

/ m ^ m7r • m7r ( cosh yiy cosh /i 2 y 

rotx = -(t^31,r)|x=f = 2^ sin 'lT — r 6 r 6 

n odd a 1 \ COsh Ml 2 COsh M2 5 



The expression in the previous parentheses is an even function of y, so it can be expanded 
in a Fourier cosine series: 



where 



x-y m7r • m7T ( 

rotx = 2^ f'lm sin — I 2L, D 

modd a * \ n odd 



12 nm COS 



nny 



(16) 



4 rf / cosh my cosh y\ nny 

D "™ = iL (^Tf ” j cos ~r dy 



Ann , Tin 
sin 



Mi M 2 



fc2 2 (m? + (tt) 2 ) + <¥) 2 ) 



The moment at the edge y = - is: 



My ~ D(—W 3 l,yy ~VW 3l , XX ) \ y _ b = ^ C 1 m D ^ + /Xjj COS 



mnx 



= ^ Flm COS 



mnx 



Therefore 



f lm — 



Fir 



D(-ci + /‘l) 



(17) 



The solution sinusoidal in the y-direction which has zero displacement on the edges and zero 

6 . 

moment M y along y = ±- is : 
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nny ( cosh V\x cosh v 2 x 



W32= L C2 " COS b V cosh V, S 



Cosh U'2 % 



w 



here v x and v 2 are roots of (14). 



Now t he rotation at the edge x — — is 



, / m ^ nny ( v\<i u 2 a 

"otx = — (n; 32 t x)|x=f = 2 ^ 2 n cos ~ ( — rq tanh — - + v 2 tanh 



= Y. c ‘» 



cos 



6 

niry 



D 



22 n 



( 18 ) 



The rotation at the edge y = - is given by: 

^ — c nir . utt / cosh t/j x cosh i/ 2 x' 

ro^y = -(i«32.v) v= | = 2 ^ C 2 n -r- sin — : 7 ; 7 

2 nodd h - V coshi/ i2 coshi/2 2 , 



where 



= E^si„^(^D 

nodd ^ ^ \m odd 



77Z7TJ* 



21 nm COS 



(19) 



D 



2 1 mn — 



4 77 Z 7 T . mn 


v \ ~ v\ 


2 SU1 O 

a 1 2 


le+fr) 2 ) e+w 2 )] 



The moment at the edge x = - is 



Mx F ( tC?32,xx ^ / ^32iyy ) |x= ^ 



nny 



V C 2n D 7q + ^ 2 ) cos 1 

^ odd ^ 

v- ~ ^Try 

2^ F 2n cos — — 



Therefore 



C 2n = 



F: 



2 n 

r>(-^r+ I'D 



( 20 ) 



In terms of the Fourier coefficients of quantities rotx , and 7 *oiy, A/ x , and A/ y , we have the 
following relationships: From equations (16), (17), (18) and (20) we have: 



11 



[rotx] = [DM 12] * [M y ] + [DA/22] * [M x \ 



( 21 ) 



u here 

( DM 12) nm = 

(DA/22) nm = 

Similarly from (15), (17), (19) and (20) we have: 

[roty] = [DA/1 1] * [M y \ + [DA/21] * [M x ] (22) 

where 

U7T . n7T 1 

DA/21 ) nm = D 2 inm * — * Sin — * — — r 77, 

6 2 D(^-/if) 

( DM 1 1 ) nm = Dum * n . 2 1 2 \ * ( Identity ) nm 

U \^ l 2 — L l l) 

Equations (21 )-(22) can be written in the matrix notation: 



' DM 11 


DM 12 ' 




' My ' 




roty 


DA/21 


DA/22 




M x 




rotx 



D 



1 2 mn 



mn . mn I 

* * sin — * — — 5 yr 

Dim * — — * ( Identity ) 



D{4 ~ "?) 



[DM] * [ M ] = rot 

Here [DM] is the flexibility matrix and the inverse [DA/] -1 is the stiffness matrix. To take 
advantage of the diagonal matrices DA/11 and DA/22, we write the following relations: 

[A/ x ] = [M] -1 * ([rotx] — [DM21] * [D Mil] -1 * [roty]) 

[M y ] = [DA/11] -1 * ([roty] — [DM12] * [M x ]) (2:1) 

[A/] = [DA/22] - [DA/21] * [DA/11] -1 * [DA/12] 
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\tiei ( omputing the vec tor solution to equations (23), we have Fourier coefficients of \l x and 
.\/ 7 and from ( 17) and (20) we can compute the Fourier coefficients of the displacement 
( or responding to this bending moment. Then uq + um — is the solution which has zero 
rotation on all the edges as well as zero displacement. 

Additional Membrane Stresses 

I ho governing equation for the additional membrane stresses in a plate with no imperfection 
is from (!)) 



Et 



-Vbf) 



“\*y * 



w 



yy 



with the boundary conditions 



(24) 



On .r — : 


fL 


d y 


= 0, 


fit 


(±i-y) dy 


— 0. u — constant 




2 






2 


(25) 


On y = : 


sh ^,xx 
2 


p, i|) dx 


- 0, 


iL *xy 


p, ±q dx 


= 0, v — constant. 



Fhe basic idea behind the solution method for calculating the in-plane Airv stress function 
( 1> is the use of double Fourier expansions for both <I> and the out of plane displacement w. 
lu order to solve the governing equation (24), we need to calculate the Fourier coefficients 
of ir from the solution in the previous subsection. The Fourier expansions for w and <1 have 
the following form: 



= E E w 

l even J even 

. = £ E . 






finx\ fjny\ 

cos v~ J ros (— J 



(26) 



l even J even 

With these expansions the boundary conditions (25) for 4> will be satisfied. The derivatives 
of tr are calculated by direct differentiation of (26), and the products of derivatives are 
evaluated at each grid point. Then the discrete double Fourier cosine transform is used to 
obtain the Fourier coefficients of the right-hand side of (24), which are stored in the Gaussian 
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curvature matrix /\ . The Fourier coefficients of <I> then follow from 



<j> = £"<» 

*0 



-j for i = 2,4 j = 2 . 4 

,(¥)' + (f)') 

The in-plane displacements u and v are obtained from the strain-displacement relations 



An appropriate plane stress solution is added to satisfy the boundary conditions on u and v 
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RESULTS 



Normal Displacement and Bending Stresses 

We lirst solved the linearized equations (7)-(S) with no imperfection and the loosely-clamped 

boundary conditions (4). The shapes of the deflected midsurfaces and values of the nondi- 

mensionalized longitudinal and transverse bending stresses for plates of various aspect ratios 

a 

and axial load ratios are shown in the figures. Note that for large - the normal displacement 

b 

and stresses on lines y = constant are nearly uniform at a distance from the ends greater 

F 

than the width of the plate. Note that as the axial load ratio — is varied from positive to 

i cr 

negative. the waviness of the plate near the ends increases. 

Maximum normal displacements, and maximum and minimum values of the longitudinal 
and transverse bending stresses are shown in separate graphs. The location of the maxi- 
mum longitudinal bending stresses is always on the centerline y — 0, while the minimum 
longitudinal bending stresses are always at the endpoints x = +— , ?/ = 0. The maximum 

transverse bending stresses arc' always on the centerline y — 0, while the minimum transverse 

b . F 

bending stresses are always on the edgeliues y = + — . Note that as the axial load ratio — 

2 l 1 ci’ 

is varied from positive to negative, the maximum stresses increase and the minimum stresses 
usually decrease. Also, note that for most cases the minimum transverse bending stress is 
the greatest stress in absolute value. 

F 

All values of — were greater than the buckling loads, given in the following table: 

r a * 



a/l. 


i 


1.2 


1.4 


1.6 


2 


4 


CO 


Huckliug 
Load F/ F cr 


-2.52 


-2,15 


-2.41 


-2.22 


-1.99 


-1.81 


-1.75 



F 

As a. check on these results, the maximum normal displacements for the case = 0 

/ cr 

(located at the origin x = y = 0) were checked against Roark and Young (1975), the 

transverse bending stress distribution on the line x = 0 was verified to be nearly parabolic 

for large -, and the buckling loads were checked against Timoshenko and Gere (1961). 
b 
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Longitudinal Bending Stress M x 6/(pb ), a/b = 1, F/F 
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Longitudinal Bending Stress M x 6/(pb ), a/b = 2, F/F, 
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Longitudinal Bending Stress M x 6/(pb ), a/b = 16, F/F, 
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Longitudinal Bending Stress M x 6/(pb ), a/b = 1, F/F ( 
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Longitudinal Bending Stress M x 6/(pb ), a/b = 4, F/F, 
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Longitudinal Bending Stress M x 6/(pb 2 ), a/b = 16, F/F 
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Longitudinal Bending Stress M x 6/(pb ), a/b = 1, F/F ( 
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Longitudinal Bending Stress M x 6/(pb 2 ), a/b = 16, F/F 
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Transverse Bending Stress M v 6/(pb ), a/b = 2, F/F, 
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Transverse Bending Stress M v 6/(pb ), a/b = 4, F/F, 
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Transverse Bending Stress M v 6/(pb ), a/b = 16, F/F 
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Transverse Bending Stress M v 6/(pb ), a/b = 1, F/F 
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Transverse Bending Stress M v 6/(pb ), a/b = 2, F/F 
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Transverse Bending Stress M y 6/(pb^), a/b = 4, F/F, 
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Transverse Bending Stress M v 6/(pb ), a/b = 16, F/F ( 
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Transverse Bending Stress M v 6/(pb ), a/b = 1, F/F 
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Transverse Bending Stress M y 6/(pb 2 ), a/b = 4, F/F 
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Transverse Bending Stress M v 6/(pb ), a/b = 16, F/F, 
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Normal Displacement wEt /(pb ), a/b = 1, F/F ( 
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Normal Displacement wEt /(pb ), a/b = 4, F/F ( 



i 

\ 



LO 



cd 

^ton 

C\J C\J i- O 
O O O 

dod 



h- co 

CD xfr 
CD CO 
O O 
O O 

d d 




CO 


LO 


C\J 


LO 


T— 


LO 


O LO 


o 


C\J 


o 


T— 


o 


o 


o 


o 


o 


o 


o 


d 


o 


o 




d 




d 




d 


d 



42 



Normal Displacement wEt /(pb ), a/b = 16, F/F, 
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Normal Displacement wEt /(pb ), a/b = 1, F/F 
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Normal Displacement wEt /(pb ), a/b = 4, F/F, 
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Normal Displacement wEt /(pb ), a/b = 16 , F/F, 
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Normal Displacement wEt /(pb ), a/b = 1, F/F ( 
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Normal Displacement wEt /(pb ), a/b = 2, F/F, 
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Normal Displacement wEt /(pb ), a/b = 4, F/F, 
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Normal Displacement wEt /(pb ), a/b = 16, F/F, 
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Maximum Longitudinal Bending Stress vs Axial Load and Aspect Ratio 
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Maximum Transverse Bending Stress vs Axial Load and Aspect Ratio 
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Minimum Transverse Bending Stress vs Axial Load and Aspect Ratio 
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Maximum Normal Displacement vs Axial Load and Aspect Ratio 
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- 0.5 



Additional Membrane Stresses 



We next studied the importance of nonlinearities by solving equation (9) with no imperfection 
and the clamped boundary conditions (4) or (5), The graphs show the nondimensionalized 
maximum and minimum additional membrane stresses due to bending. Note that the ad- 
ditional membrane stresses are at least an order of magnitude smaller than those in the 
preceding section from the linear equations. Note also that the main effect of preventing 
lateral displacement of the edges y = ±- is to increase the additional transverse membrane 
stress, whereas the additional longitudinal membrane stress is unaffected. 

For these graphs the following parameter values which are typical of ship plating were 
chosen: 

E = 30,000,000 psi 
v = .3 
b = 36 in 
L — .5 in 
p = 18 psi 

The additional membrane stresses corresponding to other parameter values may be obtained 
from these same graphs with the vertical axis scaled in accordance with (10). 
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F/F, 



Imperfections 



Finally, we solved the linearized equations (7)-(8) with no pressure, an imperfection, and the 
simply-supported boundary conditions (6). Let us consider imperfection shapes of the form 

inn x tiny 

wq = cos cos — — with m — 1,3,5,... and n — 1,3,5,... 

a b 

or 

tnnx . nny 

Wq — sin sin — - — with 'm = 2,4,6,... and u = 2,4,6,... 

a b 

where m and n are the number of half waves in the longitudinal and transverse directions, 
respectively. For these shapes we obtain the analytical solution 



w — 



Wq 



F rr ( mb n 2 a x 2 



— + 



F V 2a 2 mb. 



+ 1 



The bending stress resultants then follow from (3) hy differentiation: 



M, = 



n 2 D 



a 2 b 2 



Wq 



2 \ 2 
n a 



For I mb 
F I 2 a 2mb 



+ 1 



etc. 



F 

Of course, to prevent buckling we require — - > — 1. 

*cr 

An alternate imperfection shape proportional to the deflection of the clamped plate under 
small uniform pressure without inplane load was also allowed in the Fortran code, but for 
brevity we do not report on this case here. 
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CONCLUSIONS 



We have developed a Fortran code for calculating the displacements and stresses in rect- 
angular plates subjected to the effects of combined loading, various boundary conditions, 
and imperfections. The next step is to extend the code to a three-dimensional assembly 
of plates and beams more representative of ship grillages. Progress in the analyses of such 
structures has recently been made by Danielson tt al (1988, 1990, 1993, 1994). 



ACKNOWLEDGEMENTS 

The authors were supported by the Naval Surface Warfare Center (C'arderock Division) 
and by the Naval Postgraduate School Research Program. Dr. D. P. Kihl of the Survivability, 
Structures, and Materials Directorate at NSWC proposed the title problem and offered 
several valuable suggestions. 



67 



REFERENCES 



Bud. M. I). and Steele, C. R. (1991). Separated solution procedure for Lending of circular 
plates with circular holes. ASME Applied Mechanics Reviews 44. no. 11, part 2, S27-S35. 

Bird. M. D. and Steele, 0 R. (1992). A solution procedure for Laplace's equation on multiply 
connected circular domains. ASME Journal of Applied Mechanics 114. part 1. 398-404. 

Danielson. D. A., and Hodges, I). H. (1988). A beam theory for large global rotation, 
moderate local rotation, and small strain. ASME Journal of Applied Mechanics 55, 179-184. 

Danielson. D. A., Kihl, D. P., and Hodges, D. H.(1990). Tripping of thin-walled plating 
stiffeners in axial compression. Thin-Walled Structures 10, 121-142. 

Danielson. D. A., Cricelli, A. S., Frenzen, C. L. and Vasudevan, N. (1993). Buckling of 
stiffened plates under axial compression and lateral pressure. International Journal of Solids 
and Structures 30, 545-551. 

Danielson. D. A. ,(1994). Analytical tripping loads for stiffened plates. International Journal 
of Solids and Structures. To appear. 

Hughes, 0. (1983). Ship Structural Design. Wiley. 

Kang, L. C., Wu, C. H., and Steele, C. R.(1993). Fourier series for polygonal plate bending: 
A very large plate element. To appear. 

Kwok, F. M. W., Kang, L. C., and Steele, C. R.(1991). Mode III fracture mechanics analysis 
with Fourier series method. ASME Applied Mechanics Review 44, 166-170. 

Roark, R. .J., and Young, W. C.(1975). Formulas for Stress and Strain. McGraw-Hill. 

Sikora, J. P., Dinsenbacher, A., and Beach, J. E.(1983). A method for estimating lifetime 
loads and fatigue lives for swath and conventional monohull ships. Naval Engineers Journal 
63-85. 



68 



Szilard, R. (1974). Theory and Analysis of Plates. Prentice- Hall. 



Timoshenko. S. P.. and Gere, J. M.(1961). Theory of Elastic Stability. McGraw-Hill. 

Timoshenko, S. P., and Woinowskv-I\rieger( 1959). Theory of Plates and Shells. McGraw- 
Hill. 



69 



DISTRIBUTION LIST 



Director (2) 

Defense Tech Information Center 
Cameron Station 
Alexandria, VA 22314 

Research Office (1) 

Code 81 

Naval Postgraduate School 
Monterey, CA 93943 

Library (2) 

Code 52 

Naval Postgraduate School 
Monterey, CA 93943 

Professor Richard Franke, Code MA/Fe (1) 

Department of Mathematics 
Naval Postgraduate School 
Monterey, CA 93943 

D.A. Danielson, Code MA/Dd (30) 

Department of Mathematics 
Naval Postgraduate School 
Monterey, CA 93943 

Dr. Marc Lipman (1) 

Mathematical, Computer and Information 
Sciences Division 
Office of Naval Research 
800 Quincy Street 
Arlington, VA 22217-5000 

Dr. David Kihl (2) 

Survivability, Structures, and Materials Directorate 
Naval Surface Warfare Center-Carderock Division 
Bethesda, MD 20084-5000 



UUULtY KNOX I IRPadv 




3 2768 00331679 



